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Abstract. We discuss the possibility whether superheavy elements can be produced in Nature by the astro- 
physical rapid neutron capture process. To this end we have performed fully dynamical network r-process 
calculations assuming an environment with neutron-to-seed ratio large enough to produce superheavy nu- 
clei. Our calculations include two sets of nuclear masses and fission barriers and include all possible fission 
channels and the associated fission yield distributions. Our calculations produce superheavy nuclei with 
A ~ 300 that however decay on timescales of days. 

PACS. 26.30.Hj r-process - 27.90. -l-b Properties nuclei with A > 220 - 25.85.-w Fission reactions 



1 Introduction 

More than 40 years ago, the possible existence of an is- 
land of superheavy nuclei was proposed (see e.g. Refs. [TJ 
013]), due to the appearance of shell closures of protons 
and neutrons at mass numbers beyond the then available 
nuclear data. Various models predicted shell closures at 
Z = 114 and N — 184 3 . 4 . g).[6|[7] . More recent micro- 
scopic calculations [SlI^l fTUlfTTUT^] suggest that the magic 
proton number could be higher than Z = 114 (possibly 
Z=120, 124 or 126) and the neutron shell closure occurs 
at N = 172 or 184. In fact, shell closures are somewhat 
washed out to regions of shell stabilization around these 
numbers [13 . Shell closures make superheavy nuclei stable 
against the dominating decay channel in this region: spon- 
taneous fission. The highest fission barriers are expected 
for nuclei just above the closed proton shell (Z = 114) and 
just below the closed neutron shell at N = 184. This is of 
high relevance for the decay half-life of beta-stable nuclei 
in this mass region [T4llT51IT51IT71[I51 ■ 

Over the last thirty years, several experimental cam- 
paigns at Berkeley, GSI, RIKEN, and JINR Dubna have 
explored this predicted island of superheavy elements of 
increasingly stable nuclei around atomic number 114 (see 
e.g. [11l2CjIl21l2^l33] ;) and have found by now elements 
with charge numbers up to Z = 118 [24 l l2l p 6] . The decay 
chains observed so far suggest the existence of a region 
of enhanced stability created by shell effects. Alpha-decay 
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chains with fission, occurring at the end of the chain, ex- 
plored the boundary of the island at low Z and N. The 
question remains now how far the island extends in the 
other directions (low Z, high N; high Z, low N; high Z, 
high N) and what experimental paths can lead us to these 
regions. 

A related question has been addressed almost as long, 
i.e. whether superheavy elements are/could be produced 
in Nature [27Jf2S]. The only candidate process for astro- 
physical nucleosynthesis of these heavy elements is the 
r(apid neutron capture) process (for reviews see [19,30, 
I31|b Alternative terrestrial, non-accelerator, experiments 
have been discussed [32] . [33 i [34 l f35 | f36 | . 

While the very early discussion made it plausible that 
superheavy elements can be produced in r-process nu- 
cleosynthesis [57], more systematic approaches, including 
neutron-induced fission estimates, predicted that the r- 
process path would reach areas dominated by fission and 
the path to superheavy elements would be blocked in na- 
ture, dependent, however, on mass model uncertainties |37[ 

The first extended table of fission barrier predictions 
by Howard and Moller [Uj permitted to engage in more 
systematic calculations. Thielemann, Metzinger and Klap- 
dor 42J performed beta-delayed fission calculations for a 
number of mass models in combination with the above fis- 
sion barrier predictions. In their application to r-process 
calculations [33], they came to the conclusion that the 
r-process path ends in a 100% beta-delayed fission re- 
gion and superheavies cannot be reached. Their results 
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were based on a combination of the Hilf mass formula 44 
(which contains a slightly too steep mass parabola) and 
the Howard and Moller |41j fission barriers, which by now 
turn out to be somewhat underestimated [45] . The combi- 
nation of both effects led to an overestimation of the effect 
of fission 29 in the r-process path around Z = 92. 

Since then, after considering these deficiencies, stud- 
ies aimed at using r-process nuclei for cosmochronometry 
have neglected fission [46] and did not address the ques- 
tion of production of superheavy nuclei. In recent years 
improved predictions of fission barrier tabulations for ex- 
tended ranges of nuclei have become available [47,48,15, 
HH], being based on the Extended Thomas-Fermi (ETFSI), 
Thomas- Fermi (TF), Finite Range Droplet Model (FRDM) 
or the Hartree-Fock Bogoliubov (HFB) approach. This 
also led to new neutron-induced fission and beta-delayed 
fission calculations [STJIIFl] that permitted to address the 
formation of superheavy elements and long-lived cosmochro 
ters. In order to explore these issues we have performed a 
new set of r-process calculations with few preliminary re- 
sults reported in 52 . 53 , 54 . In this manuscript, we present 
detailed calculations aiming to answer the following ques- 
tions: (i) is the r-process experiencing strong fission effects 
during the build up of heavy nuclei by neutron captures 
and beta-decays or can unstable extremely neutron-rich 
nuclei with mass numbers up to A — 300 and beyond be 
produced? (ii) Does the subsequent decay to beta-stability 
pass by areas of the nuclear chart where fission dominates, 
hence blocking the production of long-lived superheavy 
nuclei? (iii) Can this area of fission dominance be sur- 
passed during the r-process towards higher charge num- 
bers, so that subsequent beta-decay ends with finite abun- 
dances in the valley of beta-stability, followed by alpha- 
decay chains towards the island of superheavy nuclei? This 
three options/questions are indicated in Fig. [T] 

Our paper is structured in the following way. Section 
2 discusses the nuclear input used: masses and fission bar- 
riers of nuclei far from stability. Section 3 discusses the 
r-process model and presents the results of the nucleosyn- 
thesis calculations. Section 4 analyzes the results address- 
ing the reliability of the nuclear input due to comparison 
with data. Section 5 presents the conclusions. 



2 The Role of Masses and Fission Barriers 

If one wants to address the question whether the island 
of superheavy elements was reached in natural neutron 
capture environments, like final products of r-process nu- 
cleosynthesis, a complete knowledge of all possible reac- 
tion and decay properties is required. This includes neu- 
tron capture, photo-disintegration, neutron-induced fis- 
sion, beta-decay, beta-delayed fission, spontaneous fission 
and for sufficiently large temperatures gamma-induced fis- 
sion. In order to provide this information in unknown 
territory, theoretical prediction are required. Prediction 
of these properties rely therefore on masses (to deter- 
mine reaction Q- values), fission barriers, optical potentials 
for transmission coefficients of particle channels, trans- 
mission coefficients for gamma-transitions (determined ei- 




Fig. 1. Possible pathways of the r-process in the nuclear chart 
ndiseiid line). On the path and the subsequent decay back to 
beta-stability there exist three options to encounter or avoid 
fission on the way to superheavy nuclei (options i)-iii) in the 
text), in option iii) super heavies would be made by producing 
progenitors with even higher charge numbers which then de- 
cay by subsequent beta- and alpha-decays. The present inves- 
tigation, limited to network calculation up to Z=110, cannot 
explore this latter option. 



ther from giant resonance properties or microscopic ap- 
proaches), level densities of excited states (obtained from 
back-shifted Fermi gas or combinatorial approaches) [49] 
[STl[S3[Si[SZl[SHllSSll^[SIl[ni[S3]- Spontaneous fission has 
been described with different levels of sophistication in- 
cluding: a simple formula fitted to experimental data |64[ 
[55] . macroscopic-microscopic approaches [TSIISBIIBT] . and 
Skyrme- Hartree-Fock approaches [T8"ll4"9"] . 

For a number of years extended and increasingly re- 
liable mass predictions have been available, the Finite 
Range Droplet Model (FRDM) [58], the Extended Thomas 
Fermi Model with Strutinski Integral (ETFSI) [69], the 
"microscopic" Duflo-Zuker mass formula [70] , and the Har- 
tree-Fock-Bogoliubov approach with Skyrme functionals [TO 
172] . Relativistic and non-relativistic mean field approaches 
which do not yet reach an accuracy level in ground state 
masses of 0.6-0.7 MeV have not been considered in r- 
process simulations. 

Global predictions of fission barriers have been scarce. 
After the Howard and Moller [IT] predictions (which are 
now known to produce too small barriers), new attempts 
have only been performed in the last decade [15,18,47,48, 

They are based on the ETFSI model @7J, the Thomas- 
Fermi approach [48] , with a similar macroscopic mass para- 
bola as FRDM, combined with shell corrections of the Fi- 
nite Range Liquid Drop model FRLDM [T5] - similar, but 
not identical to FRDM, and the Skyrme-Hartree-Fock- 
Bogoliubov model [49,63: based on the density functional 
used for the HFB- 14 masses [74]. The Skyrme Hartree- 
Fock approach has also been used for a sensitivity study 
of fission barriers for superheavy nuclei [T5J and for the de- 
pendence of the fission barrier on excitation energy [73] . 
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Fig. 2. Contour plots of the largest fission barrier heights from 
the TF (top) [48] and ETFSI (bottom) ,47 mass models. This 
figure also indicates to which extent data sets are available in 
the nuclear chart. HFB barriers are only publicaly available in 
the rage Z = 90-102 [49] and FRDM barriers [T5| only close 
to stability. 



Fig. [2] shows contour plots of the largest fission barrier 
for individual nuclei from the TF (top) and ETFSI (bot- 
tom) models. The calculation of beta-delayed or neutron- 
induced fission requires the knowledge of mass differences 
that should be derived consistently from the same model 
as the barriers 0B]. Unfortunately this is not always pos- 
sible, however based on the similarity of the underlying 
model we expect that the most reasonable combinations 
for barriers and masses are ETFSI/ETFSI, TF/FRDM, 
FRLDM/FRDM, and HBF-14/HFB-14. Of these four sets 
only two (ETFSI/ETFSI, TF/FRDM) can currently be 
used for global r-process calculations which require knowl- 
edge of nuclear properties from stability to the neutron 
drip line for an extended range of Z values. The Bruslib 
database of astrophysical reaction rates []] includes fission 
barriers based on the HFB- 14 model for the range Z = 90- 
110 and neutron induced rates up to Z = 102. The recent 
FRLDM barriers QH cover the range Z = 78-125 but 
always close to the beta-stability line. 



1 http : //www. astro .ulb . ac.be/bruslib 



Fig. 3. Contour plots of Bf — S n for the combinations 
TF/FRDM (top) and ETFSI/ETFSI (bottom). While there 
exist positive values of Bf — S„ for neutron numbers around 
N = 184 also for Z > 95 in the bottom figure (indicating that 
neutron-induced fission is not permitted), this is not the case 
for TF/FRDM. The neutron-drip line is shown for the appro- 
priate mass modes (FRDM and ETFSI). It can be seen that 
in the vicinity of the neutron-drip line neutron-induced fission 
is clearly inhibited. 



In figures [3] and |4] we present the different reaction and 
decay properties for neutron-induced, beta-delayed and 
spontaneous fission. Figure [3] shows contour plots of the 
difference between fission barrier and neutron separation 
energy, Bf — S n . Nuclei for which this value is below 2 MeV 
are expected to fission immediately after capturing a neu- 
tron. Given the fact that fission barriers do not suffer from 
strong odd-even effects while the neutron separation ener- 
gies are affected by those, we expect that neutron-induced 
fission will mainly occur for even-N r-process nuclei. 

Figure [4] shows the dominant decay channel (alpha- 
decay on the neutron-deficient side of stability, beta-decay, 
beta-delayed fission in regions where the beta-decay Q- 
value is larger than the fission barrier of the daughter 
nucleus, and spontaneous fission), see also 50 , 5"Tl [54l[75] . 
As the TF barriers provide only information about the 
height of the largest barrier we have estimated the spon- 
taneous fission half-lives using the simple parametrization 
of Ref. [64], but with coefficients adjusted to the avail- 



4 



I. Petermann et al.: Have Superheavy Elements been Produced in Nature? 



N 



115 
110 
105 
100 
95 
90 
85 

115 
110 
105 
100 
95 
90 
85 



' 1 1 1 I 

beta delayed fission 

spontaneous fission 

a decay 



TOO 



120 



140 



160 180 200 220 240 
N 



beta delayed fission 
spontaneous fission 
a decay 

P" -= 



00 



120 



140 



160 180 
N 



200 



220 



240 



Fig. 4. Dominant decay channels: alpha-decay on the proton- 
rich side of stability - not of importance here, beta-decay, beta- 
delayed fission in regions with positive Qp — B / (daughter) and 
regions where spontaneous fission half-lives are shorter than 
beta-decay half-lives. 



able experimental fission half-lives related to the TF fis- 
sion barriers: 



log*i/ 2 (s) = 8.08% - 24.05, 



(1) 



where the barrier Bf is given in MeV. We use the same 
parameters for the ETFSI barriers so that for the same 
barrier both models predict the same spontaneous fis- 
sion half-life. Alpha-decay half-lives are obtained using the 
Viola-Seaborg semi-empirical relationship between half- 
life and alpha decay energy with parameters determined 
in ref. [7B]. We use beta-decay half-lives and beta-delayed 
neutron emission probabilities from ref. |77| which are cur- 
rently the only values available for the mass range of in- 
terest here. As these data only extend to Z = 110, this 
is the upper limit for our r-process network. Whenever an 
experimental half-life is known we use this value in our 
calculations. 

Beta-delayed fission rates are determined based on the 
FRDM beta-decay rates [77] using an approximate strength 
distribution for each decay based on the neutron-emission 
probabilities. For each fissioning nucleus and fission chan- 
nel the fission yields, including the amount of neutrons 



produced, are determined |78j using the statistical code 

abla cigisia. 

While there are obvious differences for fission barrier 
predictions of the different models shown in Fig. [2] they 
become more prominent in Bf—S n (Fig. 3| and the promi- 
nent decay channels as shown in Fig HT especially close 
to the neutron shell closure N = 184. The content of 
Figs. [3] and [4] provides the information for judging which 
decay/reaction channel is dominant in a particular region. 
However, the determination of the neutron-induced fission 
rate requires also the knowledge of the neutron density, 
n n , of an r-process environment. One important result is 
that Bf — S n is always positive close to the neutron-drip 
line, independently of the mass model. This is expected 
as the fission barriers increase while the neutron sepa- 
ration energies decreases towards the neutron drip line. 
This means that an r-process does not experience strong 
neutron-induced fission if proceeding close to the drip line 
for N < 184. However, at the N = 184 magic shell the 
drip line moves closer to beta-stability and the quantity 
B f - S n becomes smaller. For the TF/FRDM model, its 
value is around 2 MeV for nuclei around Z ~ 94, N ~ 186 
for which the neutron-induced fission dominates over the 
(n, 7) channel. This inhibits the production of heavier nu- 
clei, due to neutron-induced fission. The ETFSI/ETFSI 
model predicts values of Bf — S n « 8 MeV in the same 
region, remaining large along the r-process path to heav- 
ier nuclei. The build up of heavy nuclei continues and 
only ends by beta-delayed fission in the region Z ~ 105, 
N ~ 220, see black region in figure [3] Notice that both 
models predict very similar neutron separation energies 
but rather different fission barriers. 

In the range of elements up to Z = 115, regions dom- 
inated by spontaneous fission will always be encountered 
during the decay to beta-stability. For the TF/FRDM case 
this is unavoidable, as an extended and connected region 
dominated by spontaneous fission has to be passed. In 
the ETFSI/ETFSI case, the spontaneous fission region is 
split into two "islands" divided by a range of about 10 
units in neutron number around N = 184, where beta- 
decay dominates without being followed by delayed fission. 
Thus, beta-decay back to stability will have the chance to 
proceed through the channel between the islands, until 
it encounters the region with N < 184. Dependent on the 
beta-decay half-lives encountered, superheavy nuclei could 
exist for an extended period of time before fissioning. We 
cannot make statements for a possible production of nuclei 
beyond Z — 110 as this is the upper limit of our present 
network calculations due to the restricted availability of 
the nuclear input. For this region extended sets for theo- 
retical predictions have been published [181 after finishing 
our calculations and in the future we plan to explore their 
impact in r-process nucleosynthesis. 

The decay and reaction inputs discussed in this sec- 
tion and displayed in Figs. [2j [3] and [4] provide the basis 
for applications to the r-process calculations given in the 
following section. 
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3 r-Process Calculations with fission 

3.1 Astrophysical Sites 

A high availability of neutrons, leading to a large ratio of 
neutrons to seed nuclei (typically Fe-group nuclei or be- 
yond, up to A — 80) can occur in two types of explosive 
astrophysical environments: (a) Very neutron-rich mate- 
rial, which is obtained under high nuclear matter/neutron 
star densities, where electron capture on protons by degen- 
erate electrons with high Fermi energies cause a proton- 
to-nucleon ratio Y e of about 0.1, is ejected and leads to an 
r-process during the expansion phase. Possible examples 
are the ejecta from neutron star collision events, [81,82. 
83,84], or other environments ejecting highly neutronized 
material like fast rotating core collapse supernovae with 
strong magnetic fields and jet ejecta [81 ^|8"6"1I5T] ■ (b) Mat- 
ter which is only slightly neutron-rich but experiencing a 
fast expansion and high entropies. In such a case the reac- 
tions 4 He(cva, n) 9 Be(a, n) 12 C and/or 4 He(aa, 7) 12 C are 
responsible to move matter from 4 He to 12 C followed by 
a fast sequence of reactions producing heavier seed nuclei 
and few neutrons remaining. At high entropies, oc T 3 /p, 
three body reactions are suppressed due to the low density 
and/or high number of photons hindering the production 
of heavy seed nuclei and result in large amounts of 4 He 
(with N — Z) and a tiny amount of heavy seed nuclei. 
This results in large remaining neutron-to-seed ratios. Ex- 
amples of this case are the neutrino-driven wind [88,89, 
I901I91U92|[9^11M1193] from the nascent proto-neutron star 
after supernova core collapse (provided that it can pro- 
duce neutron-rich ejecta [961197] ) and matter ejected from 
accretion disks around black holes [98 , 99 . 

In the above two options, the high neutron-to-seed ra- 
tio is the dominant driving force behind a rapid neutron 
capture (r-)process. However, dependent on the environ- 
ment, the densities and temperatures can vary and cause 
differences in the exact working of the r-process. In or- 
der to produce nuclei as heavy as Uranium and Thorium 
which are only produced by the r-process, one requires 
conditions with neutron-to-seed ratios of the order of 150. 



3.2 r-Process model calculations 

The previous subsection outlined the still existing uncer- 
tainties in astrophysical sites/conditions which will then 
also be affected further by nuclear uncertainties discussed 
in section [2] We stress that, independent of the actual en- 
vironment, the decisive feature of an r-process to produce 
superheavy elements is a sufficiently high neutron-to-seed 
ratio. In the following, we assume assumed conditions cor- 
responding to fast expanding ejecta with high entropy, in 
order to ensure a large enough neutron-to-seed ratio. 

We start our calculations at a temperature of 10 GK 
for which Nuclear Statistical Equilibrium is applicable and 
assume that the matter follows a homologous adiabatic 
expansion where the density behaves as: 

p(t)=poexp(-t/T). (2) 



We use t = 3 ms following the hydrodynamical simula- 
tions of ref. [100] . At later times hydrodynamical simula- 
tions show that the evolution is not anymore homologous 
and can be approximated by 101 -,1021 1103] : 

pe>-»(^)'. m 

where the parameter A represents the time scale on which 
the matter evolves from conditions of almost constant den- 
sity (t « A) to constant velocity (t » A) 95J. We use 
a value of A = 2 s in agreement with hydrodymanical 
simulations 100, 35] • 

In order to attain a large neutron-to-seed ratio we 
assume a constant entropy of 200 fc/nucleon and initial 
Y e — 0.35 that results in a neutron-to-seed ratio of 290 at 
3 GK. The temperature is determined by the condition of 
constant entropy using the equation of state of ref. [104] . 
Our network code is based on the XNet code of ref. [105] 
that has been extended to treat implicitly all fission reac- 
tions. 




Time (s) 



Fig. 5. Evolution of temperature and density for a hot and 
cold r-process depending on the choice of ti in Eq.(3). Notice 
that for both cases the n/seed ratio drops below 1 in the range 
t=l-2s. 

As discussed in section 3.1, the r-process can occur 
under quite different astrophysical conditions. In the fol- 
lowing, we present two extreme cases for the late time evo- 
lution (see figure [5]) by choosing t\ as the time for which 
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the temperature has reached values of 1.15 and 0.1 GK, 
respectively. In the following, we refer to these two cases 
as "hot" and "cold" r-process. 

The nucleosynthesis in these two cases differs in two 
major aspects. First, in the "hot" r-process calculations, 
we find that there is a continuous production of seed nu- 
clei during the whole duration of the r-process as charged- 
particle reactions never completely freeze-out (see ref [106 ). 
This is not the case in the "cold" r-process, as for such 
low temperatures charged-particle reactions are too slow 
when compared with the dynamical timescale of r-process 
nucleosynthesis. The continuous production of seed nu- 
clei manifests itself in enhanced abundances for nuclei in 
the range A — 20-90. This is demonstrated in figure [6] 
that shows the final r-process abundances based on the 
TF/FRDM input. Calculations with the ETFSI/ETFSI 
input show the same features. 




90 120 150 
Mass Number A 



240 



Fig. 6. Final abundances for the hot and cold r-process sce- 
narios after 1 Gy; i.e. after the decay of unstable isotopes (see 
text for details). 



The second difference relates to the role of photodis- 
integration reactions during the r-process. At the larger 
temperatures of the hot r-process the photodisintegration 
reactions are fast enough to establish a (n, 7) (7,^) 
equilibrium. This chemical equilibrium between both re- 
action types along an isotopic chain, defines a narrow "r- 
process path" , with a nucleus of dominant abundance in 
each chain, acting as waiting point for beta-decays. This is 
the "classical" r-process picture as introduced in ref. |1071 

una. 

In the "cold" r-process, photodisintegration is too slow 
to establish an equilibrium with neutron captures. The 
evolution is rather determined by neutron captures and 
competing beta-decays which opposite to photo-dissociation 
do not vary much from isotope to isotope under r-process 
conditions. As a consequence several nuclei in each iso- 
topic chain are substantially populated which leads to a 
broader r-process path. Calculations in the cold r-process 
scenario have been already presented in 109,I 110| I95|. how- 
ever, with different aims than discussed here. 

Besides this more general discussion, addressing the 
working of the r-process, we want to concentrate now on 



the focus of the present investigation, the effect of fission 
and the path to superheavy nuclei. For this purpose we 
have performed calculations in the hot and cold r-process 
scenarios using the TF/FRDM and ETFSI/ETFSI input 
sets. We focus in the following on the mass regime A > 180 
that is relevant for our discussion. The abundance dis- 
tribution at the moment where the neutron-to-seed ra- 
tio reaches 1, marking the beginning of neutron-capture 
freeze-out, is shown in figs. [7] and [8] In each case a com- 
parison is given between a hot and a cold r-process, shown 
in Fig. [7 for a nuclear input set TF/FRDM for fission bar- 
riers and masses and in Fig. [§] for the set ETFSI/ETFSI. 

We see that in both cases (hot and cold r-process) and 
for both sets of nuclear input the abundances decline sub- 
stantially for nuclei A > 280, extending in the cold case for 
a few more mass units. As discussed by 38,109 in a hot 
r-process (n, 7) <=^ (7, n) equilibrium ensures a substantial 
leakage into the fission channel, hindering the production 
of heavier elements once an isotopic chain is reached where 
neutron-induced fission can occur. This occurs in the re- 
gion Z = 95-100 close to or slightly beyond N = 184, 
where as seen in Fig. [3] Bf — S n < 2 McV. Once neutrons 
are exhausted, Fig. B] (top) shows that spontaneous fission 
can become the dominating fission channel. In the cold r- 
process scenario the evolution is dominated by individual 
nuclei with low fission rates, allowing the production of 
heavier nuclei than in the hot scenario. 

To explore the dependence on the nuclear input, figs. [9] 
and [lO] show the evolution of fission rates for the different 
channels for a hot and cold r-process using the fission rates 
and masses TF/FRDM and ETFSI/ETFSI. There are sev- 
eral general features that are worth to be mentioned. Dur- 
ing the r-process phase, defined by a neutron-to-sccd ratio 
larger than one, neutron-induced fission dominates. Once 
the neutron-to-seed becomes smaller than one and the r- 
process freeze-out starts, the rate of neutron-induced fis- 
sion suddenly decreases. At this moment, the r-process 
material starts to beta-decay to stability, producing beta- 
delayed neutrons via beta-delayed neutron emission that 
can induce neutron-captures but also neutron-induced fis- 
sion. The latter is particularly important for the matter 
accumulated at JV = 184 (A ~ 280) that during the 
beta-decay to the stability feeds the region with Z ~ 95, 
N ~ 175 where neutron-induced fission can operate again, 
see fig. [3j This revival of neutron-induced fission is seen in 
all the calculations by a second hump that is delayed by 
the time scale of the successive beta decays. Once neutron- 
induced is revived it sustains itself by a mechanism sim- 
ilar to a chain reaction and continues to be the domi- 
nating channel till the neutron density becomes too low 
at times around several tens of seconds when beta-delayed 
fission and spontaneous fission dominate. Spontaneous fis- 
sion, which is dominant for nuclei closer to stability, is al- 
ways less important during the working of the r-process, 
but becomes larger in the freeze-out phase, when the path 
moves closer to stability. In the cold case the path is more 
spread out in abundances, as discussed above, and en- 
counters, therefore, the spontaneous fission region (closer 
to stability) earlier, already before the decay back to sta- 
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Fig. 7. Results from calculations for hot (top) and cold (bottom) r-process conditions, utilizing fission barriers and mass 
predictions from the models ETF and FRDM. In both cases the abundances are shown as a function of A as well as in terms of 
a contour plot, where the abundances of nuclei are indicated by their color. The abundances are given at the point of neutron 
freeze-out, i.e. when the ratio of neutrons to heavy nuclei has dropped down to 1. That means that on average each nucleus will 
not capture more than one neutron after this point in time. 
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fission barrier/mass model selection TF/FRDM. 
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the different channels shown for a hot (top) and cold (bottom) 
r-process, but with the fission barrier/mass model selection 
ETFSI/ETFSI. 



bility sets in. For nuclei with the highest mass numbers 
(see also Fig.Q beta-delayed fission can also be impor- 
tant at late times. 

While for TF/FRDM nuclear input we have seen that 



lower of the two spontaneous fission islands, at a charge 
number Z > 110, i.e. outside the limits of the presently 
used nuclear network. 



Fig. [TT] shows (for the ETFSI/ETFSI combination of 
neutron-induced fission prevents the build-up of (super) heavyfission barriers and masses) results of hot and cold r- 
elements, the situation looks different when applying ETFSI/process calculations at that point in time when Z=109 



ETFSI for fission barriers and masses (see Fig. |8[). This 
is understandable from Fig. |3j where we see that the re- 
gion with Bf — S n < is further removed from the drip- 
line (and r-process path) than for the TF/FRDM case. As 
beta-delayed fission is always less prominent than neutron- 
induced fission during the r-process build-up, we also ex- 
pect matter to continue to proceed to elements heavier 
than those included in the present calculation (see sect. [2]). 
The relative role of the three fission channels shown in 
Fig. [lO] is similar to that shown in Fig. [9j but with the 
exception that neutron-induced and beta-delayed fission 
do not fully hinder the production of nuclei heavier than 
those included in the present network, limited to Z < 110. 
Another difference is that spontaneous fission becomes 
more prominent in the late freeze-out phase when the 
path moves closer to stability. This can be understood 
as neutron-induced and beta-delayed fission do not fully 
prevent the production of heavy nuclei. 

It was already realized from Fig. [4] that the area where 
spontaneous fission dominates is split for the ETFSI/ETFSI 
case, leaving a "free passage" in the vicinity of N = 184 
[52"] . This permits beta-decay in the direction of the valley 
of stability. The decay path will encounter, however, the 



nuclei are still present or just decayed to Z=110, i.e. when 
reaching the limits of the nuclear network in the late decay 
phase. In both calculations this happens after about 12h. 
Thus, when utilizing ETFSI/ETFSI fission barriers and 
masses, such (superheavy) nuclei exist for this amount of 
time. From Fig. [4] we know, however, that they will finally 
also reach the lower island of spontaneous fission and not 
reach the valley of stability as final destination. 

As noticed above in the ETFSI/ETFSI case the r- 
process involves nuclei heavier than those included in our 
network and consequently we cannot determine the an- 
swer to question (iii) of the introduction, i.e. the possibil- 
ity of producing nuclei heavy enough to circumvent the 
region of fission dominance also during the decay to sta- 
bility. This option is not available for the TF/FRDM case 
as sufficiently heavy nuclei are never produced during the 
r-process. 



4 Discussion and Summary 

In this paper we examined whether superheavy elements 
can be synthesized in nature by the r-process. As the r- 
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Fig. 11. Top and bottom figures show the abundance patterns of hot and cold r-process calculations with ETFST/ETFSI 
fission barriers and masses. The plots refer to the latest point in time when still Z=109 nuclei are present, or just have decayed 
to Z—100. These are the results of decay via the "free passage" region of spontaneous fission shown in Fig. [4] 



process conditions in astrophysical explosions are not fully 
understood, we utilized two sets of calculations, relating 
to the so-called hot and cold environments, where either 
photo-disintegrations play or do not play an important 
role. For conditions with sufficient neutron-to-seed ratios 
to produce actinide and heavier nuclei, we came to essen- 
tially two types of conclusions. 

For one set of fission barriers and mass predictions 
(TF/FRDM) the r-process is terminated by neuton- induced 
and beta-delayed fission before superheavy nuclei can be 
reached. For the other set (ETFSI/ETFSI) neutron-induced 
and beta-delayed fission do not prevent the build-up of su- 
perheavy nuclei even extending beyond Z > 110. Unfortu- 
nately, this is the edge of our network calculations, related 
to the current limits of available global nuclear input. In- 
dependent of this shortcoming, our calculations predict 
that nuclei with masses of the order A — 290 are produced 
and survive for about 12 hours after an r-process event 



(which only takes seconds). After several beta decays also 
these nuclei will encounter regions of spontaneous fission 
an decay to medium mass nuclei before reaching the valley 
of stability. 

The present investigation is another example showing 
the importance of nuclear input for r-process nucleosyn- 
thesis. A reduction of the underlying uncertainties by ei- 
ther experiment or theory is desirable. This should also 
include an extension of the calculations to nuclei with 
higher charge, Z > 110, and neutron number to deter- 
mine if the r-process can produce progenitors with such 
high mass numbers that during the decay to the stability it 
can circumvent the region of fission dominance predicted 
by current models (N ~ 184, Z > 100). A first step to- 
wards this goal has been achieved in a recent study |18j . 
which shows that within the uncertainties among differ- 
ent Skyrme functionals it is likely that spontaneous-fission 
will hinder the production of elements beyond Z = 120. 
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